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Finite-Length Multi-Input Multi-Output 
Channel Shortening Pre-filters 

Related Application: 

5 This application claims priority from Provisional application No. 60/158,713 

filed on October 8, 1999. This application is also related to a Provisional 
application No. 60/158,714, also filed on Octobers, 1999. 

Background of the Invention: 

10 The combination of maximum likelihood sequence estimation (MLSE) with 

receiver diversity is an effective technique for achieving high performance over 
noisy, frequency-selective, fading channels impaired by co-channel interference. 
With the addition of transmitter diversity, the resulting multi-input multi-output 
3 (MIMO) frequency-selective channel has a significantly higher capacity than its 
ms single-input multi-output (SIMO) or single-input single-output (SISO) counterparts. 
^; The use of maximum likelihood multi-user detection techniques on these 
3 frequency-selective MIMO channels significantly outperforms single-user detection 

techniques that treat signals from other users as colored noise. However, MLSE 
™J complexity increases exponentially with the number of inputs (or transmit 
IMo antennas) and with the memory of the MIMO channel, making its implementation 
% over sever inter-symbol interference (IS!) channels very costly. 
3 The MIMO channel can be modeled as a collection of FIR filters (i.e., an 

FIR filter between each input point (e.g., transmitting antenna) and each receiving 
point (e.g. receiving antenna), and the "memory of the channel" corresponds to the 
25 number of taps in the FIR filters. 

The Discrete Matrix Multitone (DMMT) was shown to be a practical 
transceiver structure that asymptotically achieves the MIMO channel capacity 
when combined with powerful codes. It uses the Discrete Fourier Transform 
(DFT) to partition the frequency responses of the underlying frequency-selective 
30 channels of the MIMO systems into a large number of parallel, independent, and 
(approximately) memoryless frequency subchannels. To eliminate inter-block and 
intra-block interference, a cyclic prefix whose length is equal to the MIMO channel 
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memory is inserted in every block. On severe-ISI MIMO channels, the cyclic prefix 
overhead reduces the achievable DMMT throughput significantly, unless a large 
FFT size is used which, in turn, increases the computational complexity, 
processing delay, and memory requirements in the receiver. 
5 In short, the computation complexity increases exponentially with the 

number of taps in the FIR filters that may be used to model the channel. 

N. Al-Dhahir and J. M. Cioffi, in "Efficiently-Computed Reduced-Parameter 
Input-Aided MMSE Equalizers for ML Detection: A Unified Approach," IEEE Trans. 
Information Theory, pp. 903-915, May 1996, disclose use of a time-domain pre- 
10 filter in the receiver to shorten the effective channel memory and hence reduce the 
cyclic prefix overhead and/or the number of MLSE states. The disclosed 
approach, however, is for SISO systems, and not for MIMO systems. 

Summary 

45 An advance in the art is realized with a MIMO pre-filter that is constructed 

f; from FIR filters with coefficients that are computed based on environment 
3 parameters that are designe^chosen. Given a transmission channel that is 

modeled as a set of FIR filters with memory v , a matrix W is computed for a pre- 
3 filter that results in an effective transmission channel B with memory Nb, where 
[20 A/^ < V , where B is optimized so that B^^ = argmin^ trace{R^^) subject to selected 

3 constraints; Ree being the error autocorrelation function. The coefficients of W, 
which are sensitive to a variety of designer constraints, are computed by a 
processor within pre-filter at the front end of a receiver and loaded into an array of 
FIR filters that form the pre-filter. 

25 

Brief Description of the Drawing 

FIG. 1 shows the major elements of a receiver in accord with the principles 
disclosed herein; 

FIG. 2 presents the structure of pre-filter 30; and 
30 FIG. 3 is a flowchart describing the method carried out by processor 220 

within pre-filter 30. 
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Detailed Description 

FIG. 1 depicts the general case of an arrangement with riq transmitting 
antennas 11-1, 11-2, ... Il-ziq, that output signals (e.g., space-time encoded 
5 signals) to a transmission channel, and /?o receiving antennas 21-1, 21-2, ... 21-a?o. 
Each transmitting antenna p outputs a complex-valued signal Xp, the signals of the 
A?q antennas pass through a noisy transmission channel, and the rio receiving 
antennas capture the signals that passed through the transmission channel. The 
received signals are oversampled by a factor of / in element 20 and applied to pre- 
10 filter 30. Thus, the sampling clock at the output of element 20 is of period Ts=T/l, 
where 7 is the inter-symbol period at the transmitting antennas. The transmission 
channel's characterization is also referenced to Tg. In the illustrative embodiment 
i disclosed herein, therefore, pre-filter 20 develops n\ output signals that are applied 
fi to a conventional multi-input receiver 40, and the received signal can be 
expressed by 

3 yi^^=i^;/^■^'xi^.+"i^^ (1) 

3 where y^/^ is the signal at time k at the receiving antenna, /?*|'^^ is the m^^ 
i coefficient (tap) in the channel impulse response between the transmitting 
U antenna and the/^ receiving antenna, and n^^^ is the noise vector at the/^ 
yo receiving antenna. The memory of this path (i.e., the largest value of ah for which 
/)^'^^ is not zero) is denoted by v^''^\ It not unreasonable to assume, however, that 
the memory of the transmission channel is the same for all ij links (n, x such 

links), in which case v*''^^ = v . Alternatively, the v^''^^ limit in equation (1) can be 
set to that v which corresponds to maximum length of all of the xn^ channel 

25 input responses, i.e., v = max^^ v^''^^ . All of these variables in equation (1) are 

actually /x1 column vectors, corresponding to the / time samples per symbol in 
the oversampled FIG. 1 arrangement. By grouping the received samples from all 
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10 



15 



20 



A7o antennas at symbol time k into an nJx^ column vector y^^ , one can relate y^^ 
to the corresponding n, x1 (column) vector of input samples as follows 

yk = T.K^R-m+n,, (2) 

m=0 

where Hm is the MIMO channel coefficients matrix of size njxrif , x^_^ is a size 
n^x^ input (colunnn) vector, and is a size n^/xl vector. 

Over a block of symbol periods, equation (2) can be expressed in matrix 
notation as follows: 
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or, more compactly, 



)fk+Nf-\k ~ ^^k+Nf-tk'V ^k+Nf-tk ■ 



(4) 



The subscripts in equation (4) indicate a range. For example k-\-Nf-^ :k 
indicates the range k + Nf-^ and k , inclusive. 

It is useful to define the following correlation matrices: 

^x/ ^ ^[^k+Nf-^k'vyk+Nf-^:k] ~ ^xx'^ (^) 
Ryy = ^[y/f+A/,-1:/fy/f4-A/,-t/f ] = ^R^^H + R^^ , (6) 
^xx = ^[^k+Nf-tk-v^k+Nf-l-k-v] C^) 
^nn = E[^k+Nf~tk^k+Nr-tk ] ■ 

It is assumed that these correlation matrices do not change significantly in time or, 
at least, do not change significantly over a time interval that corresponds to a 
TDMA burst (assumed to be much shorter than the channel coherence time), 
which is much longer than the length of the pre-filter, in symbol periods denoted by 
A/f. Accordingly, a re-computation of the above matrices, and the other parameters 
disclosed herein, leading to the computation of pre-filter coefficients, need not take 
place more often than once every TDMA burst. 
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Once H, Rxx and Rnn and known, Rxy and Ryy are computed by R^^H* and 
HR^^H* + R^^ , respectively. 

Given the MIMO channel matrix H with + 1 members (Ho,H^,"-HJ, the 
objective is to create a MIMO pre-filterW (element 30 in FIG. 1) with Nf matrix 
5 taps, i.e., matrix W = [w^ ^Nr-^J ' that equalizes H so as to create an 

overall transmission channel for receiver 40 that corresponds to a matrix B with 
memory N^^ , where N^«v . 

The matrix B can be expressed as B = [Bq B^ • • ■ B^^^ ]^ where each Bi 
is of size x . 

10 The MIMO channel-shortening pre-filter W (element 30) is conditioned, or 

:J adjusted, to minimize the equalization Mean Squared Error (MSE), defined by 
n MSE = /face(Rgg) , where Re© is the autocorrelation matrix of the error vector Ek 
0 that is given by 

3 ^k^Nf-tk~v ~ ^ Vk+Nf-U ' (9) 

15 where the augmented MIMO matrix, B* , is 

:i B*-[o,,,,^, b; b; ... b;^ o,^,,^,]-[o,^,,^, b* o^,,^^], (io) 

y A is the decision delay that lies in the range 0<{N^ -\-v-N^-T), and 
3 s = Nf+v-Nf^-A-"l, The rij x tij error autocorrelation function Ree can be 
expressed by the following: 

Ree -^l^k^k] 

=B (R^^-R^^R;;jR^)B ^^^^ 
= B*R^B 
= B*RB, 

where R is a sub-matrix of R-^ determined by A. 

Using the orthogonality principle, which states that E[E^y*^;y^_^.;,] = 0 it can 

be shown that the optimum channel-shortening pre-filter and target impulse 
response filters (W and B, respectively) are related by 
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**opf ~ '^opt^xy^yy 

= B^pfR^^H (HR^^H -^Rnn)'^ (12) 

The last line shows explicitly that the MIMO channel-shortening pre-filter consists 
of a noise whitener R;;^ , a MIMO matched filter H*, and a bank of FIR channel- 
shortening pre-filter elements. 

5 It remains to optimize B such that the MSE is minimized, which may be 

obtained by computing the parameters of B that, responsive to specified 
conditions, minimizes the trace (or determinant) of Ree- The following discloses 
two approaches to such optimization. 

Under one optimization approach the coefficients of B are constrained so 
f:3iO that some coefficient of B is equal to the identity matrix, I. A solution subject to 
f S this Identity Tap Constraint (ITC) can be expressed by 
I B^J = argmine trdce{RJ subject to - 1,^ , (13) 

j; where <D* ^ [o,^,,,, 1,^ 0.,..,(A/,-m)] and 0 < m < A/, . It can be shown that the 

optimum MIMO target impulse response and the corresponding error 
autocorrelation matrix are given by 
ry Bj^ = R-^<I>(c|)*R-^c|))-^ and (14) 

I RS.in=(0*R-^<I>)-\ (15) 

As indicated above, R is affected by the delay parameter A . Unless dictated by 

the designer, the delay parameter A , which can range between 0 and 
20 (Nf + v-N^-^), is chosen to minimize the trace of R^J^^nm • Similarly, the index 

parameter m, which ranges between 0 and , and which that affects O , is 

chosen to minimize the trace of R JJ^,-^ . 

Under a second optimization approach the imposed constraint is B*B = 1^ . 

A solution subject to this Ortho-Normality Constraint (ONC) can be expressed by 
25 B^p^^ = argming frace(R^J subject to B*B = l„ , (16) 

Defining the eigen-decomposition 
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R - UZU* = Ud/ag(c7o.aT-<T,^(yv,.iH)U* , (17) 
where (Jq > cr^ • • • > c7^^(a/,^.ih , then the optimum MIMO target response and the 
resulting error autocorrelation matrix are given by 

^Opt "* ^/7,(A/^+1)-l] ' (18) 

5 where e/ is unit vector with a 1 at position /, and O's elsewhere, and 

Re^min = ^'^S^(c^/.,A/, ^ ' % Cr^,(A/,+1H) • (19) 

Illustratively, if a?/=3 and Nb=3, B^^^^ =U[eQ,e^Q,e,,], meaning that B^^^^ is a three 

column matrix comprising the 9^^ through the lf ^ columns of matrix U. Stated in 
words, the optimum MIMO target impulse response matrix is given by the n\ 
10 eigenvectors of R that correspond to its n\ smallest eigenvalues. The delay 
parameter A (0<A<A/^4-v-A/^-1) that affects R is optimized to minimize the 
trace (or determinant) of R^j„ . 

With the above analysis in mind, a design of a prefilter 30 can proceed for 
: any given set of system parameters, which includes: 

l5 • MIMO channel memory between the input points and the output point of the 
actual transmission channel, v , 
The number of pre-filter taps chosen, A/f, 
The shortened MIMO memory, A/b, 
The number of inputs to the transmission channel, riu 
20 • The number of output derived from the transmission channel, /lo, 
The autocorrelation matrix of the inputs, Rxx, 
The autocorrelation matrix of the noise, Rnn, 
The oversampling used, /, and 
The decision delay, A . 
25 The structure of pre-filter 30 is shown in FIG. 2. In the illustrated 

embodiment, it comprises two main components: processor 220 and filter section 
210. 

Filter section 210 in the FIG. 2 illustrative example comprises a collection of 
FIR filters that connect the rio input array of signals from sampling circuit 20 to an 
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H/ output array of points. That is, there are jxi FIR filters P/,/ , that couple input 

point j to output point /. 

Processor 220 is responsive to the no signals received by antennas 21 and 

sampled by circuit 20, and it computes the coefficients of W, as disclosed above. 
5 Wo is a matrix that defines the coefficients in the 0^^ tap of the j x / FIR filters, Wi 

is a matrix that defines the coefficients in the 1®* tap of the / x / FIR filters, etc. 
The method of developing the parameters of pre-filter 30, carried out in 

processor 220, is shown in FIG. 3. Block 100 develops an estimate of the MIMO 

channel between the input points and the output point of the actual transmission 
10 channel. This is accomplished in a conventional manner through the use of 

training sequences. The estimate of the MIMO channel can be chosen to be 

limited to a given memory length, v , or can be allowed to include as much 
3 memory as necessary to reach a selected estimate error level. That, in turn, 

depends on the environment and is basically equal to the delay spread divided by 
m 7s. 

2 Following step 100, step 110 determines the matrices, Rnn, Rxx, Rxy, and 

Ryy. The matrix Rnn is computed by first computing n = y-Hx and then computing 

^ the expected value E[n*n] see equation (8) above. The matrix Rxx is computed 

from the known training sequences -- see equation (7) above - (or is pre- 
fo computed and installed in processor 220). In may be noted that for uncorrelated 
^ inputs, Rxx=l. The matrices Rxy and Ryy are computed from the known training 
sequences and the received signal or directly from H and Rnn- see equations (5) 
and (6) above. 

Following step 110, step 120 computes R"^ =R;,;, -R^^R;;jR^ , and the sub- 

25 matrix R . From equation (10) is can be seen that R is obtained by dropping the 
first 111 A rows and the last n^s rows of R-^. 

In accordance with the ITC approach, selecting some value of 0 < m < A/^ 
allows completion of the design process. Accordingly, following step 120, step 
130 chooses a value for m, develops <P* ^ [o,^,,,, 1,^ ^n^-rr^)] ^^^^''^^ ^ut 
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the computation of equation (13). Step 140 finally develops the coefficients of 
matrix W in accordance with equation (12), and installs the developed coefficients 
within filter 210. 

In accordance with the ONC approach, step 130 computes the matrix U in a 
5 conventional manner, identifies the unit vectors e,, and thus obtains the matrix B. 
Step As with the ITC approach, step 140 develops the coefficients of matrix W in 
accordance with equation (12), and installs the developed coefficients within filter 
210. 

It should be understood that a number of aspects of the above disclosure, 
10 for example, those related to the ITC constraint and to the ONC constraint, are 
merely illustrative, and that persons skilled in the art may make various 
modifications that, nevertheless, are within the spirit and scope of this invention. 
!{ For example, the pre-filter described above generates a multi-output signal, with 
fi the number of outputs being n,, that being also the number of transmitting 
45 antennas 1 1 . This, however, is not a limitation of the principles disclosed herein. 
* The number of pre-filter outputs can, for example, be larger than n,, for example as 
3 high as /7,(A/6 + 1) . The performance of the receiver will be better with more filter 
3 outputs, but more outputs require more FIR filters, more FIR filter coefficients, and 
fj correspondingly, a greater processing power requirement placed on processor 
m 220. 
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Claims: 

1= A receiver responsive to an rio plurality of antennas comprising: 

a pre-filter having an n^xn^ plurality of FIR filters, each responsive to a 

signal that is derived from one of said rio antennas and applied to an input point, 
and each developing an output signal that contributes to one of /?/ pre-filter 
outputs; and 

decision logic responsive to said tii outputs. 

2. The receiver of claim 1 further comprising a sampling circuit interposed 
between said Hq plurality of antennas and said pre-filter that samples received 

signal at rate = j ^ where I is an integer and T is symbol rate of a transmitter 

whose signals said receiver receives. 

3. The receiver of claim 2 where />1 . 

4. The receiver of claim 1 where coefficients of said FIR filters are 
computed in a processor in response to a block of Nf symbols. 

5- The receiver of claim 4 where said processor is part of said pre-filter. 

6. The receiver of claim 4 where said coefficients of said FIR filters are 
computed once every time interval during which transfer characteristics of said 
transmission channel, H, are substantially constant. 

7. The receiver of claim 6 where said processor installs computed 
coefficients of said FIR filters in said FIR filters following each computation. 
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8. The receiver of claim 1 where said FIR filters form an array of filters, and 
said array has one FIR filter connected between each of said no input points and 
said A7/ outputs. 

9. The receiver of claim 8 where said no plurality of antennas receive 
signals, via said transmission channel, from a transmitter having a multiple number 
of transmitting antennas. 

10. The receiver of claim 9 where said transmitter has n, transmitting 
antennas. 

11 . The receiver of claim 10 wherein said decision logic is adapted to 
receive from said transmitted signals that were encoded in a space-time encoding 
schema. 

12. The receiver of claim 2 where said plurality of FIR filters is expressed 
by matrix W, and W is computed by W*^^ = B'^^R^^R^ , 

=b;,,R,,H*(HR..H*+R„J-\ or W;, =B;,(R;:+H*R-;:HrH*R-i, where R,, is 

an autocorrelation matrix of a block of signals transmitted by a plurality of 
transmitting antennas to said Ho antennas via a channel having a transfer 
characteristic H, R„„ is an autocorrelation matrix of noise received by said plurality 
of no antennas during said block of signals transmitted by said transmitting 
antennas, R^ =RxxH*, Ryy =HR^H + R„^ , and 

B*pj is a sub-matrix of matrix B'^j, where B^^, = argming trace{R^) subject to a 
selected constraint, Ree being the error autocorrelation function. 

13. The receiver of claim 12 wherein said plurality of FIR filters are 
subjected to designer constraints relative to any one or a number of members of 
the following set: transmission channel memory, size of said block, effective 
memory of the combination consisting of said transmission channel and said pre- 
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filter; n-,, no, autocorrelation matrix Rxx, autocorrelation matrix Rnn, value of factor / 
in said sampling circuit, and decision delay. 

14. The receiver where said matrix W is expressible by 

5 \N = [ Wq W, • • . W^^_i , where matrix W^, is a matrix that specifies q**" tap 
coefficients of said FIR filters. 

15. The receiver of claim 12 where said constraint restricts B so that 
B*q) = I , where <D* ^ [0„,,„,^ Kn.iN^-m)] and m is a selected constant. 
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16. The receiver of claim 15 where B = R"^<P(<P*R"^<P)"^ , R is a sub-matrix 



of a matrix R^ = R^ -R^^R^^R^ . 



17. The receiver of claim 12 where said constraint restrict B so that 



:15 B B = L . 



5 18. The receiver of claim 17 where B = U[e„;^ ••• e„_(^^i)_^],each 

y element ep is a vector having a O'element in all rows other than row p, at which 

3 row the element is 1 , and U is a matrix that satisfies the equation R = UIU* , Z 

20 being a diagonal matrix. 
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Abstract 

A multi-input, multi-output pre-filter improves operation of a multi-input 
receiver by shortening the effective memory of the channel with a set of FIR filters. 
The coefficients of these FIR filters can be fashioned to provide a variety of 
controls by the designer, for example, the value of the effective memory. 
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IN THE UNITED STATES 
PATENT AND TRADEMARK OFFICE 

Declaration and Power of Attorney 

As a below named inventor, I hereby declare that: 

My residence, post office address and citizenship are as stated below next to my name. 

I believe I am an original, first and sole inventor of the subject matter which is claimed 
and for which a patent is sought on the invention entitled Finite-Length Multi-Input Multi- 
Output Channel Shortening Pre -filters the specification of which was filed on 10/8/99, as 
application Serial No. 60/158,713. 

I hereby state that I have reviewed and understand the contents of the above- identified 
specification, including the claims, as amended by an amendment, if any, specifically referred to 
in this oath or declaration. 

I acknowledge the duty to disclose all information known to me, which is material to 
patentabiUty as defined in Title 37, Code of Federal Regulations, 1.56. 

I hereby claim foreign priority benefits under Tifle 35, United States Code, 119 of any 
foreign application(s) for patent or inventors' certificate listed below and have also identified 
below any foreign application for patent or inventors' certificate having a filing date before that 
of the application on which priority is claimed: 



I hereby claim the benefit under Title 35, United States Code, 120 of any United States 
application(s) listed below and, insofar as the subject matter of each of the claims of Uns 
application is not disclosed in the prior United States application in the manner provided by the 
first paragraph of Title 35, United States Code, 112, we acknowledge the duty to disclose all 
information known to us to be material to patentability as defined in Titie 37, Code of Federal 
Regulations, 1.56 which became available between the filing date of the prior application and the 
national or PCT international filing date of this application: 

US Provisional application 60/158,713, filed October 8, 1999. 

I hereby declare that all statements made herein of my own knowledge are true and that 
all statements made on information and belief are believed to be true; and finther that these 
statements were made with the knowledge that willfiil false statements and the like so made are 
punishable by fine or imprisonment, or both, under Section 1001 of Title 18 of the United States 
Code and that such willfid false statements may jeopardize the vaKdity of the appUcation or any 
patent issued thereon. 
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I hereby appoint flie following attomey(s) with full power of substitution and revocation, 
to prosecute said application, to make alterations and amendments therein, to receive the patent' 
and to transact all business in the Patent and Trademark Office connected therewith: 



Samuel H. Dworetsky 


(Reg. No. 27873) 


Thomas A. Restaiio 


(Reg. No. 33444) 


Michele L. Conover 


(Reg. No. 34962) 


Cedric G. DeLaCruz 


(Reg. No. 36498) 


Rohini K. Garg 


(Reg. No. 45272) 


Thomas M. Isaacson 


(Reg. No. 44166) 


Benjamin S. Lee 


(Reg. No. 42787) 


Robert B. Levy 


(Reg. No. 28234) 


Susan E. McGahan 


(Reg. No. 35948) 


Gary H. Monka 


(Reg. No. 35290) 


Jeffrey M. Navon 


(Reg. No. 32711) 


Stephen K. Pentlicki 


(Reg. No. 40125) 


Alfred G. Steimnetz 


(Reg. No. 22971) 







I also appoint the followuig as associate attomey(s), with full power to prosecute said 
plication, to make alternations and amendments therein, and to transact aU business in the 
Patent and Trademark Office connected therewith: 

Henry T. Brendzel (Reg. No. 26,844 
William Ryan (Reg. No. 26,844) 

Please address all correspondence to Henry T. Brendzel, P.O. Box 574, Springfield, NJ 
07081 . Telephone calls should be made to Hemy T. Brendzel at (973) 467-2025. 



Full name of sole inventca": Naofel Al-Dhahir 

Inventor's signature AI ' ^' " A^<DJ\^ Date l"? 1^^ 

Residence: Monistown, Morris County, NJ 
Citizenship: Iraq 

Post Office Address: 1 05 Harter Road 
Morristown, NJ 07960 



